Theory of the Pseudospin Resonance in Semiconductor Bilayers 
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The pseudospin degree of freedom in a semiconductor bilayer gives rise to a collective mode 
analogous to the ferromagnetic resonance mode of a ferromagnet. We present a theory of the 
dependence of the energy and the damping of this mode on layer separation d. Based on these 
results, we discuss the possibility of realizing transport-current driven pseudospin-transfer oscillators 
in semiconductors. 
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Introduction — The layer degree of freedom in semicon- 
ductor bilayers is often regarded [lj as an effective spin- 
1/2 pseudospin degree-of- freedom in which electrons in 
the top layer are assigned one pseudospin state, and elec- 
trons in the other layer the opposite one. In the quantum 
Hall regime, electron bilayers are sometimes [2| pseu- 
dospin ferromagnets. The appearance of these broken- 
symmetry states motivates an interest in phenomena 
which are pseudospin analogs of the very robust magneto- 
electric effects which underpin spintronics in ferromag- 
netic metals, and therefore might underpin a useful 
pseudospintronic technology. Unfortunately it appears 
likely [3| that pseudospin- ferromagnetism in semiconduc- 
tor bilayers is a phenomenon that is limited to strong 
magnetic fields or, possibly, systems with extremely low 
densities. This Letter is motivated by the observation 
that one important spintronic device, the spin-transfer 
oscillator [J], requires only collective spin-dynamics and 
not spontaneous magnetic order. In a spin-transfer oscil- 
lator transport, currents drive ferromagnetic- resonance 
collective spin-dynamics in the presence of applied fields 
strong enough to oppose hysteretic switching. Since 
pseudospin paramagnets do support a pseudospin res- 
onance collective mode, as we explain below, pseudospin 
polarized transport currents in a semiconductor bilayer 
(easily realizable using individual-layer contacting tech- 
niques Q) could, provided that the resonance is suffi- 
ciently sharp, drive collective pseudospin dynamics and 
yield a device with similar functionality. In this Letter 
we report on a theory of the damping of the pseudospin 
resonance which suggests that it is possible to design bi- 
layers with sharp pseudospin resonances. 

In the absence of tunneling a semiconductor bilayer 
supports two types of collective excitations d, Q: (i) 
an optical mode (the ordinary plasmon) with a long- 
wavelength dispersion relation oc q 1 ^ 2 and weak damp- 
ing (<x q 2 ), in which electrons in the two layers oscillate 
in phase, and (ii) an acoustic plasmon with linear dis- 
persion and strong Landau damping (oc q) by particle- 
hole excitations, in which the electrons in the two layers 



oscillate out of phase. The pseudospin resonance is a 
q = collective mode which develops from the out-of- 
phase plasmon when interwell tunneling is enabled. In 
the pseudospin language interwell tunneling favors sym- 
metric bilayer states and therefore acts like a pseudospin- 
magnetic field, which we take to act in the z direction. 
The pseudospin resonance then involves collective preces- 
sion around this pseudospin field, with y-direction pseu- 
dospins representing current flowing between the layers 
and i-direction pseudospins representing charge accumu- 
lation in one of the layers. Theoretical treatments of the 
pseudospin resonance have so far relied on the random 
phase approximation 0] (RPA), sometimes with Hub- 
bard or local-density-approximation corrections [H, 0]. 
These papers demonstrate that interactions shift the res- 
onance away from the position of the q = particle-hole 
excitations, eliminating the Landau damping process. In 
this Letter we present a theory of the pseudospin trans- 
fer resonance that is based on a systematic expansion 
in powers of the difference V— between intra- and inter- 
layer electron-electron interaction. We obtain an expres- 
sion for the interaction-induced resonance position shift 
which is exact to leading order in V- , and an expression 
for the leading order damping contributions which ap- 
pears at second order in V- . Damping of the pseudospin 
resonance is similar to damping of the ferromagnetic res- 
onance [Tot [Txj] in a metal, except that it is intrinsic and 
driven by electron-electron interactions rather than dis- 
order. The physical mechanism of damping is the pro- 
duction of two (or more) electron-hole pairs with zero 
total momentum. The phase space for these processes 
implies that the damping rate is proportional to the cube 
of the resonance frequency, implying that the resonance 
is sharp whenever its energy is small compared to the 
Fermi energy of the bilayer system. 

The model — In a semiconductor bilayer, electrons in the 
same layer interact through the two-dimensional (2D) 
Coulomb interaction V s (q) — 2-Ke 2 /(eq) (e is the dielec- 
tric constant), while electrons in different layers are cou- 
pled through the interlayer Coulomb interaction Vd(q) = 



2 



V s (q)e~ qd . We assume a spatially constant inter-layer 
tunneling amplitude which we denote by Asas/2 and 
present our theory using a pseudospin representation in 
which the tunneling term is diagonal, i.e. the represen- 
tation in which |f) refers to the symmetric combination 
of single-layer states and \\) to the antisymmetric com- 
bination. The total Hamiltonian is then (H = 1) 
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Here a is the real-spin label, a is the pseudospin label, S 
is the sample area, p q = J2k, a ,<j £ k-q/2,a, CT 5 k+q/2,Q : <T and 
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Sq—O, and, finally, 



density and the pseudospin operators (r a being Pauli 
matrices with a = x,y,z), S® ot 
V±(q) = [V s (q)±V d (q)}/2. ' 
Theory — The theory we develop in this Letter is based 
on the observation that the difference between the intra- 
and inter-layer interaction V-(q) = 7re 2 (l — e~ qd )/(eq) 
is always smaller than ire 2 d/e, which becomes a small 
perturbation when d <C max(r s aB, ab/^). Here r s — 
(7rna|) -1 / 2 is the Wigner-Seitz density parameter and 
ae = e/(me 2 ) is the Bohr radius. The above inequal- 
ity guarantees that the third term in the Hamiltonian 
(fTJ) is a small perturbation either compared to the ki- 
netic energy [~ e 2 / (r 2 eae)] which dominates in the high- 
density limit, or compared to the interaction energy 
[~ e 2 /(r s e<ZB)] which dominates in the low-density limit. 
We will therefore perform a systematic expansion for 



the pseudospin resonance frequency and damping rate 
in powers of V- (q). Our approach will be asymptotically 
exact in the limit tlCaB, and is expected to be qualita- 
tively correct for d ~ ob- 

We determine the properties of the pseudospin res- 
onance by evaluating the transverse pseudospin re- 
sponse function XS X S X (<Zi w ) — ((Sq > SZ„}}oj/S, where 
we have introduced the Kubo product {{A;B))u = 
-ilim e ^ 0+ r +oo ^e^e- £ *(* GS |[i(t), J B(0)]|* GS ) g^. 
Since our aim is to calculate the transverse mode at 
q = we will focus on the response function xs ,x S x ( tJ ) = 
XS"S x (q — 0,w) [HI]. The in-plane pseudospin operators 
satisfy the Heisenberg equations of motion, 
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S^ otl which measures the difference between charges in 
the two layers, is a good quantum number when A§as — ► 
whereas S^ ot is not conserved even in this limit because 
of the pseudospin-dependent interactions. When d — > 
these equations reduce to a pseudospin version of Lar- 
mor's theorem, in which the precession is undamped and 
its frequency is given exactly by the non-interacting par- 
ticle value Asas- 

Our theory starts by making repeated use of Eqs. © 
in the Kubo product identity [12j, [l4|: {{A;B)) U — 
(*Gs|[A-B]|*Gs)/w + i((9 t i;-B» w /w. After some alge- 
braic manipulations we arrive at the following exact ex- 
pression for X5 x 5 x ( w ) 
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is the ground-state pseudospin magnetization per unit 
area, /(k) 
fl (k) = {^ GS \S£S y 
£(k,k',w 



*GS>, 

, and 



Notice that /(k) is purely real, <?(k) is purely imaginary, 
and £(k, k',w) has both a real and an imaginary part. 
The symmetric interaction V+ does not appear explicitly 
in Eq. When V_ is set to zero (d — > 0), the inter- 
action part of the Hamiltonian is pseudospin invariant. 
Larmor's theorem then applies to the pseudospin degree- 
of-freedom and only the first term on the right hand-side 
of Eq. ([3]) survives. We refer to the Hamiltonian TL at 



V_ = as the reference system (RS), on which the per- 
turbative scheme outlined below is based. 

The key idea now is to expand xs x S x ( w ) hi powers of 
VI. For example, the ground-state pseudospin magneti- 
zation M z is expanded as M z = Mq + Ml + M% + 
where the n-th term M z n is 0(V™). The quantities /, g, 
and C are similarly expanded. Note that the zero-th or- 
der of /(k), denoted by /o(k), is a non-zero difference 
between longitudinal and transverse pseudospin struc- 
ture factors. On the other hand, the zero-th order of 
<7(k), denoted by go (k) vanishes because the RS Hamil- 
tonian is invariant under rotations by 90 degrees about 
the z-axis in pseudospin space which map S x — > S y and 

Sy ► S X . 
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The pseudospin resonance frequency is the solution of 
the equation ?R,e[xg^gx (<jj±)] = 0. The inverse of xs x s x 
can be expanded with the help of the formula X~s*S x = 

x 1 («) - xo 2 Mxi H + xo 2 M \xo 1 Hxi M - x 2 («)] + 

where Xn( w ) is the n-th order term in the expansion 
for xs x S x - To appreciate the power of Eq. we first 
use it to find u>± to first order in d = d/ciB- Keeping only 
terms up to first order in Eq. ([3|) and making use of the 
formula above for the power series of Xs^s* ( w ) we obtain 



Xs*s* ( w ) 
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A4§A SA s 
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which implies immediately that 

1 

S 2 



ALs+ 4As ^ (0) ^E/o(k) + ^). (5) 



This equation is exact to all orders in the intralayer 
Coulomb interaction V B . In the high-density (non- 
interacting) limit one can find simple analytical expres- 
sions for Mq and /o(k), Mq — (ng — n^£)j2 and 
5" 2 E k /o(k) = {Ml) 2 /2. Here n a = k 2 J(2ir) are the 
band occupation factors, k-p a being the Fermi wavenum- 
ber for band a. In this limit Eq. ([5]) simplifies to 
u)\ = A^ AS + 2A SAS XgI/_(0) + O(J 2 ). The second term, 
which supplies the interaction induced shift in the pseu- 
dospin resonance position, is a factor of two smaller than 



in RPA theory [8j. The source of this difference is easy 
to understand: our calculation includes the first-order 
exchange corrections to the resonance frequency which 
are absent in the RPA. Since V_ is independent of q at 
first order in d, corresponding to a J-function interaction 
in real space, the like-real-spin contribution to the res- 
onance position shift present in the RPA is canceled by 
exchange interactions. 

The main object of this work is to estimate the reso- 
nance decay rate, which appears first at second-order in 
V- and is zero in the RPA. (Additional interaction correc- 
tions to the resonance position 5Re(w^), which we do not 
discuss at length, also appear at second order [l5j.) The 
linewidth of the pseudospin resonance [= — 23m,(tj^)] is 
given, up to second order in d, by 



4I/ 2 (0)A SA s 
Ml 



lim 
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where £o{^) is the wavevector sum of the four-spin cor- 
relation function £(k, k', lu). This quantity can be evalu- 
ated analytically in the high-density kinetic-energy domi- 
nated limit in which it is dominated by a decay process in 
which two particle- hole pairs are excited out of the Fermi 
sea, one involving a pseudospin-flip. The second particle- 
hole excitation is diagonal in pseudospin and absorbs the 
momentum emitted by the first. We find that 



3?m £ (uj) = --(Ml) 3 6(u> - A SA s) - r^g ^ ^ S(oj - £l Q (k, k'))n k », Q n k »_ k+k /^(l - n k » +k ',/3)(l - ^k"-k,a) (7) 



where n k , a = 9(fcF Q — |k|) and fi Q (k, k') = k • k'/m + 
oAsas ■ The first term on the right-hand-side of Eq. ([7|) 
does not contribute to the life-time of the shifted res- 
onance. In Fig. Q] we illustrate the dependence of 
3m £0(10) on w - The w 3 dependence at small to is the 
double-particle-hole excitation manifestation of the fa- 
miliar Pauli-blocking reduction in the excitation density- 
of-states in a Fermi sea which underlies Fermi liquid 
theory; damping drops much more rapidly at low en- 
ergies than for ferromagnetic resonance 11| dominated 
by single-particle decay processes. Eqs. ((SJ, ([5]), ((U, and 
(|7|) constitute the most important results of this work and 
provide, to best of our knowledge, the first microscopic 
theory of the pseudospin resonance linewidth. 

Numerical results and discussion — Typical numerical re- 
sults for r_L, calculated from Eqs. © and ([7]) are shown in 
Figs. [Hand [3] In Fig.[5]we show T± as a function of Asas 



for a bilayer with density n — 8.3 x 10 10 cm~ 2 and inter- 
layer distance d = L + w = 50 A. Here L = 40 A is the 
width of each quantum well and w = 10 A is the barrier 
width (we have chosen material parameters correspond- 
ing to a GaAs/AlGaAs bilayer). The non-analytic behav- 
ior of Tj_ for Asas ~ 3 meV is due to the transition from 
the situation in which both symmetric and antisymmetric 
bands are occupied to that in which only the symmetric 
band is occupied. In Fig. [3] we illustrate the dependence 
of T± on density for a fixed value of Asas = 1-48 meV. 
Since the resonance frequency is close to Asas > these cal- 
culations predict that the pseudospin resonance can be 
very sharp, especially when Asas is small compared to 
the Fermi energy of the bilayer. On physical grounds 
we expect that the main effect of going to higher order 
in d will be to replace the bare interlayer interaction in 
Eq. ([6|) by a weaker screened interaction, further reducing 
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FIG. 1: (Color online) Imaginary part of the dynamical re- 
sponse function £o(ui) (in units of eV -1 nm -6 ) as a function 
of uj for a bilayer electron gas with n = 8.3 x 10 10 cm -2 and 
Asas = 1-48 meV. The 5-function contribution at uj = Asas 
[see first term in Eq. ([T}] has been subtracted. The red solid 



line is the asymptotic result Qm £q(uj — > oo) = 



732. 



Inset: a zoom of the low-energy region. The solid red curve is 
the expression Sim £o(u>) = — juj 3 with 7 ~ 3.41 x 10 -4 [l6| ]. 
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FIG. 2: (Color online) Intrinsic linewidth F± of the pseu- 
dospin resonance as a function of Asas f° r a bilayer with 
density n = 8.3 x 10 10 cm~ 2 and d = 50 A. The S2D curve was 
evaluated using the bare 2D interactions V B (q) and Vd{q) de- 
fined above whereas the Q2D result was evaluated with more 
realistic interactions weakened by form factors [l7| ] which ac- 
count for typical quantum well widths. 



the damping. 

Our theory of the resonance amounts to the deriva- 
tion of an anisotropic, linearized pseudospin Landau- 
Liftshitz-Slonczewski equation: 



d t M a 



d t M v = - 



AsasA^ - 
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where A4 a is the average macroscopic pseudospin po- 
larization, which becomes equal to (5 t a ot ) in the limit 
I — ► 0. In the first line of Eq. (JSj> we have added a Slon- 
czewski fl8j ] pseudospin transfer term proportional to the 
tunnel current /, which is injected in one layer and ex- 
tracted from the other. As in the ferromagnetic case, it 
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FIG. 3: (Color online) Intrinsic linewidth Fj_ of the pseu- 
dospin resonance as a function of n for a bilayer with tunnel- 
ing gap A s as = 1-48 meV and d = 50 A. The labels S2D and 
Q2D have the same meaning as in Fig. [2] 



is the reaction counterpart of the torque which acts on 
the transport quasiparticles to enable their transfer be- 
tween layers upon moving through the sample, and must 
be present because of the nearly exact conservation of 
pseudospin by interactions. In the second line of Eq. @ 
we have added a Gilbert-like damping term oc d t M. x (the 
anisotropy of the Gilbert damping in the present prob- 
lem derives from the strongly anisotropic character of the 
interaction part of the Hamiltonian) . These equations 
[which describe a damped pseudospin precession of fre- 
quency wj_ and damping rate T± about the steady state 
values M y (t -> 00) = J/(eA SAS ), M x {t -> oc) = 0] 
are similar to those which describe spin-transfer torque 
oscillators 0] in ferromagnets and suggest that simi- 
lar, and possibly more flexible, devices could be real- 
ized in semiconductor bilayers. We anticipate that the 
pseudospin resonance in ferromagnets will have negative 
rather than positive dispersion, because of the q depen- 
dence of V- (q) . The roles of this property, and the fact 
that the single-particle and collective excitation frequen- 
cies are not widely separated, are difficult to fully antici- 
pate. Nevertheless, this work suggests that experimental 
studies of non-linear transport in bilayers have great po- 
tential. 
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